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Abstract 

Animal locomotion employs different periodic patterns known as animal gaits. In 1993 Collins 
and Stewart achieved the characterization in quadrupeds and bipeds by using permutation sym- 
metries groups which impose constrains in the locomotion centre called Central Generator Pattern 
(CGP) in the animal brain. They modelled the CGP by coupling four non linear oscillators and 
with the only change in the coupling it is possible to reproduce all the gaits. In this work we 
propose to use coupled chaotic oscillators synchronized with the Pyragas method not only to char- 
acterize the CGP symmetries but also evaluate the time serie behaviour when the foot is in contact 
with the ground for futures robotic applications. 



PACS numbers: 05.45.-a 



I. INTRODUCTION 



Lately the analysis of animal gaits have regain interest in between scientist of different 
areas. The moderm analysis represents the gait as cyclic patterns in the movements of 
symmetrically placed limbs. A cycle means the interval between successive footstrikes of the 
same foot during the dynamical process. The factor of footstrike f3 is the fraction of cycle 
when the foot is in contact with the surface. For simplicity we assume (3 to be the same 
for all animal feet. The relative phase of foot (</>) is defined as the fraction of cycle between 
the contact of the surface of the foot of reference and another foot when this one enters in 
contact with surface^. In this study the relative phase plays a crucial role to formulate the 
symmetries while this is not same the factor of footstrike f3, which will no be taken into 
account in this work. The mammal phenotypes have evolved in two kind of gaits, bipedal 
gaits where the limbs can be out of phase (walking or running) or in phase (jumping or 
hopping). Quadrupedal gaits with a more complex behaviour of the realtive phase. The 
natural gaits are^: Walk, the limbs move with a quarter of cycle out of phase, where there 
is a quarter cycle phase difference in between both fore limbs as well as in between hind 
lims, and half a cycle between diagonals. Trot, the diagonal legs move in phase and this 
pair is half a cycle out of phase of the other one. Pace or Rack, the fore and hind left 
or roght limbs are paired and moved half cycle out of phase in between the pairs. Canter, 
the right front leg and the left hind leg move in phase, the left front leg and the right hind 
leg move half a cycle out of phase in between themselves and out of phase with respect to 
the first pair (it was found in horse that pattern changes from walk to trot, to canter to 
gallop, when the speed increase). Bound, the fore legs move in phase, as well as the hind 
legs, while they move half a cycle out of pahse. Transverse Gallop, the front and hind left 
(right) legs move one quarter of cycle out of phase, the front (hind) limbs are slightly out of 
phase between themselves. Rotatory Gallop, similar to the tranverse gallop but the left 
and right limbs have their pattern exchange. Pronk, the four limbs move together an in 
phase. 

Biological model assume that the animal nrevous system contains a variety of Central 
Pattern Generator^ (CPG), each oriented to specific action. For instance, the locomotion 
CPG controls the rhythm of mammal gait4, in the case of quadrupedal mammals this is 
modeled by a system of coupled cell where erach cell is composed by a set of neurons 
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directly responsible to harmonize the movement of the leg^. A simplifield mathematical 
model of locomotion CPG consists in replacing each cell by a non linear oscilator— . This 
model has been studied using differnt method: bifurcation theory, numerical simulations 
and phase repons o 7 i 8 i 9 i 10 . 

The idea to study rhythmic patterns in animal gait using symmetries models was intro- 
duced by Hildebrand^, Schoner et clIf^. The concept of symmetries in coupled cells as a 
model for locomotion CPG in quadrupedal mammals was firt used by Collins and Stewart^. 
A model for locomotion CPG for quadrupedal mammals consits in a ring of four coupled 
nonlinear oscillators. Each oscillator represents a limb of the animal. The stability and 
breakdown of the symmetries play an effective role in the validity of the model. Golubitsky 
et a££. argued that symmetries present in above model for walk, trot and pace are not 
adequade for quadrupeds, sience the trot and pace correspnd to conjugate solutions which 
have same stability and they depend on initial conditions. Many quadrupeds move with 
pace but do not trot (camels) or viceversa (horse), unless they are trained. In this work we 
propose a different coupling mechanism in order to avoid the problem of multiple conjugated 
solutions^. 

The CPG is modeled by the following system of ordinary differential equations: 



where j = 1 . . . 4 mod 4 is the index of cell, X e R n is the state vector and / : R™ — > R n 
is a nonlinear velocity vector field. We defined the symmetry of ring to the permutations of 
cell that preserves the coupling, that is to say a permutation a of {1 . . .4} numbers on the 
phase space X = (Xi, X 2 , X 3 , X±) is: 



f{X j ) + h ] {X^X ]+l ) 



(1) 



dt 



aX = (X (T -i(i),X tT -i( 2 ),X lT -i(3),X tr -i(4)) 



(2) 



then a is a symetry of ring in 



F(aX) = aF{X). 



(3) 



where F(.) = /(.) + h(.). Then we deduced the coupling conditions they are: 



hj(aXj_i,aXj + i) — h a (j)(Xj-i, Xj + i). 



(4) 
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If we defined [i,j] as the action to exchange X\ for Xj the symmetries of the ring are: 
{ [1,4]; [2,3]; [1,2]; [4,2] }. From eq-Hwe deduced hx(.) = h 3 (.) and h 2 {.) = h 4 (.). These 
symmetries are called Type-2 by Collins and Stewart^. Another kind of symetry is called 
symmetry of phase changed. Assuming that X(t) is a periodic solution with minimal period 
(cycle) T, and 7 represents the symmetry (ij) to permute Xi for Xj, then jX(t) will be 
a periodic solution if the trajectories {X(t)}t and {^X(t)} t coincide. Therefore the only 
solution is the existence of pahse slip 9 such that jXj(t) = Xj(t + 9). The pair (7,6?) is a 
spatio-temporal symmetry where 9 is a phase slip. Finally we define as primary gai1>^, those 
gaits modeled by identical output signal of each cell but out of phase. 

We associated the index of cell to each limb as follow, j = 1 hind left, j : = 2 fore left, 
j — 3 fore right and j ' = 4 hind right. The possible symmetries of the primary gait for four 
legged animals characterized by Type-2 arrays are: 

Table 1: Simmetries asocied with gait 



Gait 


Symmetry 


Group 


Stopped 


(1,9) [a, 6) ((3,9) (a(3,9) 


D 2 x S 1 


Pronk 


(1,0) (a,0) (AO) (a/3,0) 


D 2 


Pace 


(1,0) (a,\) ((3,\) (a/3,0) 


f)D 
u 2 


Bound 


(7,0) (a,0) (13,1) (a/3,i) 


Di 


Trote 


(1,0) (a,i) ((3,0) (a(3,\) 


Dk 


Rotatory Gallop 


(J.0) ((3,\) 


Z 2 


Transverse Gallop 


(1,0) (a/3, \) 


7 F 


Canter 


(1,0) 


I 



TABLE I: Where a = (12)(34), (3 = (13)(24), a/3 = (14)(23) and S 1 refers to all cyclic group of 
pahse slip mod 1. D represents the diedral subgroup and Z all the cyclic subgroups. The tilde 
indicates the existence of phase slip symmetry. The notation 1 represents a half cycle out of phase. 

Here we study the possibility to use a ring of coupled chaotic oscillators to produce the 
locomotion of quadrupeds based on the symmetries of primary gaits using Pyragas control 
theory. 
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II. A CPG MODEL 



Here we represented each cell by Rosier oscillators (see eq.-EJ) coupled using Pyragas^ 
method with random initial cnditions. We use Rossler chaotic oscillator since this is the 
only one that have shown to synchronize to simulate the primary gait. This behaviour does 
not happen for Van der Pool^ or Showaltep^ oscilators even if the single oscillator reproduce 
the necessary output. Primary gate is important for any digital aplication on mechanical 
limb. 



f 

'dt 



dy 
dt 

dz 



X 



0.2 y (5) 



§ = 0.2 + z(x-c) 

We use a direct synchronization mechanism where the master variable is "y" and other 
one are the slave varibles. We use a delay time serie for obtain delay feedback value. Then 
coupling functions are: 

hi(Xi, X i+1 , Xi-i) = h (y l - 1 (t - t) - y^t)) 

(6) 

+9i {Vi+i{t - t) -y l {t)) 
The symmetry conditions associated to a Type-2 array limits the range of ki and g\ values. 
In this case g\ = k2, #2 = #3 = k^, 9a = k\, k\ = k2 and k^ = k%. An esquematic form 
it is depected in figure-^ The delay time r G N, and the nonlinear constan c (see eq.-EJ) 
play an important role in the wave pattern obtained. We consider as output of each cell 
(oscilator) the value of the variable Xi(t), where it will be compouse by threshold function: 

, si x > 2.0 

Q(x) = { (7) 
1 si x < 2.0 

This defined a mapping from phase space into binary matrix space of 2x2. We asso- 
ciate the value "1" to the state ll limb on ground 1 and the value "0" to the state ll limb in 
movement 1 , not on the ground. Finally the matrix representation is from now on: 

Fore Left Fore Right 
Hind Left Hind Right 

Therefore the gait is nothing but the sequence of matrices of succesive states representing 
the symmetry of CPG. For instance the pronk is give by the sequence of matrices: 



C 




This allows a clearer visualization of the symmetries of primary gait. Although we lose 
the time interval between different patterns. This is highly important for application in 
robotics, not only it conditions actuator answers, also for the fact that mechanical inner 
may introduce undesirable instabilities. This is the reason that we have to analyse the time 
each pattern stays in periodic sequence. 

We consider two type of ad hoc combinations for the coupling constants, which are the 
most represntatives in between the values tried. We call SA model when k\ = k% = 0.1 and 
fc 2 = &4 = 0.001, SB a k\ = k 3 = 0.1 and k 2 = k^ = —0.001. Those values were selected 
under the assumption of strong coordination between the limbs associated to each cerebral 
hemisphere, while they are weakly correlated when they belong to different hemisphere.^ 




FIG. 1: Cell coupling diagram. 
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III. NUMERICAL RESULTS 



There is a dependence of the output, on the time delay r, which is quiet robust under 
variations of the parameter c (see eq-0), and it is independent of the model used, SA or 
SB. There is a change of state as a function of the time delay, which follow the order Chaos 
— > Periodic Oscillations — > Stable Fixed Point — > Primary Hopf Bifurcation — > Secondary 
Hopf Bifurcations to Chaos when r increases from zero. This can seen in Fig-El where we 
plot the maximum Lyapunov exponet of the systems versus r. Since it is necessary to have 
a stable time interval in between patterns, we restrict ourselves to the interval 6 < r < 38. 




FIG. 2: Maximum Lyapunov exponent verus delay. 



A. SA Coupling 

For 6 < r < 13, the periodic gait obtained is: 




which corresponds to the pronk and the symmetry is D 2 . On the other hand for r = 14 
the limit cycle is not stable any longer ant it appears an asymptotic stationary state, which 
produces the single sequence correspond to the symmetry D 2 x S 1 , therefore a stop, since 
all limbs are on ground. For r = 34, the fixed point loses its stability and it becomes 
unstable, where the orbits converge to a single limit cycle. The pattern found for the delay 
34 < r < 38 is : 






This has a symmetry D 2 , which corresponds to the gait trot. In this case each step 
involves the movement of all the limbs, gait which is only observed in havy quadruopeds, 
above 1 ton, such as girafes and buffalos^. We could not find any other patterns in this 
interval of time delay. 



B. SB Coupling 

For 6 < r < 13 the coupled chaotic system oscillate in a stable limit cycle, which in this 
case produce the patterns: 






which has symmetry D 2 and corresponds to a gait bound. We should mention that the 
true bound as the one observed in Siberian squirrel^ implies the state of all limbs on air. On 
the other hand in this case all four limbs are on the ground in order to generate another 
step which does not exist in nature. This drawback, we resolve applying "nof operator 
each matrix. As in the SA case, for 14 < r < 32 the coupled system has stable fixed point 
and the limit cycle becomes unstable, which correspond to stop. For r = 33 the fix point 
become unstable and a stable cycle appears producing the patterns: 

1 1 
1 1 

which corresponds to gait pace with symmetry D 2 . As describes above the horse never 
sets all four limbs on the ground, then applying u nof operator each matrix, we have resolved 

8 






it. For this coupling the gait is not structurally stable for all time delay values. From r = 34 
until 4 r = 37 the gait changes to another periodic pattern: 




With symmetry {(I, 0) (a, |) (/?, |)}. This symmetry does not correspond to any primary 
rait observed in nature. For r = 38 the system generates again a gaite like pace. 



IV. CONCLUSIONS 

Pyraga's direct synchronization (see eqEJ) is a novel coupling mechanism between cell for 
Type-2 networks when it is used as CPG model. Within this we can avoid the conjugate 
undesirable solutions. But no natural patterns appear in SB model. Also we can not look 
for canter and transverse gallop gaits. This network type is not adequately for a natural 
GCP model. However it is useful as an artificial CPG mechanism in robotic science. 
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